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ABSTRACT 

We present the largest, publicly available, sample of Damped Lyman-a systems (DLAs) along 
Swift discovered Gamma-ray Bursts (GRB) line of sights in order to investigate the environmental 
properties of long GRB hosts in the 2 = 1.8 — 6 redshift range. Compared with the most recent 
quasar DLAs sample (QSO-DLA), our analysis shows that GRB-DLAs probe a more metal 
enriched environment at 2 > 3, up to [X/H] ~ —0.5. In the z = 2 — 3 redshift range, despite 
the large number of lower limits, there are hints that the two populations may be more similar 
(only at 90% significance level) than at higher redshifts. Also, at high- 2 , the GRB-DLA average 
metallicity seems to decline at a shallower rate than the QSO-DLAs: GRB-DLA hosts may be 
polluted with metals at least as far as ~ 2kpc from the GRB explosion site, probably due to 
previous star-formation episodes and/or supernovae explosions. This shallow metallicity trend, 
extended now up to 2 ^ 5, confirms previous results that GRB hosts are star-forming and have, on 
average, higher metallicity than the general QSO-DLA population. Finally, our host metallicity 
measurements are broadly consistent with the predictions derived from the hypothesis of two 
channels of GRB progenitors, one of which is mildly affected by a metallicity bias, although more 
data are needed to constrain the models at 2 > 4. 

Subject headings: gamma-ray: burst - techniques: spectroscopic - quasars: absorption lines - galaxies: 
general - galaxies: ISM 


1. Introduction 

One of the fundamental aspects of the forma- 
tion of the first stars and galaxies is the actual 
conversion of the primordial hydrogen clouds into 
the first massive, almost metal free, objects (Pop- 


^NASA Postdoctoral Program Fellow, Goddard Space 
Flight Center, Greenbelt, MD 20771, USA 

^Institute for Computational Cosmology, Department of 
Physics, Durham University, South Road, Durham, DHl 
3LE, UK 

^Carnegie Observatories, 813 Santa Barbara Street, 
Pasadena, CA 91101, USA 

^Department of Astronomy and Astrophysics, 
UCO/Lick Observatory, University of California, 1156 
High Street, Santa Cruz, CA 95064, USA 

®Kavli Institute for Cosmology and Institute of As- 
tronomy, University of Cambridge, Madingley Road, Cam- 
bridge, CB3 OHA, UK 


ulation III star, Barkana & Loeb 2001). This first 
generation of stars, at 2 > 10, disappeared quite 
rapidly due primarily to strong negative feedback 
effects (Karlsson et al. 2013; Bromm 2013; Yoshida 
et al. 2008; Greif et al. 2012; Chen et al. 2014). 
However, some of these objects and the following 
generation of stars probably ended their lives in 
very energetic explosions, either as pair-instability 
supernovae or as long Gamma-ray Bursts (GRB, 
see for a review Meszaros 2013), which can be de- 
tected by current and future high-energy missions 
up to the highest redshifts. 

Thanks to the Swift satellite (Gehrels et al. 
2004), hundreds of GRBs have been discovered, 
even up to 2 ~ 8 (Tanvir et al. 2009; Salvaterra 
et al. 2009; Cucchiara et al. 2011b). These high -2 
GRBs can be used to test cosmic star- format ion 
rate models as well as the cosmological chemical 
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enrichment (Kistler et al. 2009; Robertson & Ellis 
2012; Tanvir et al. 2012; Salvaterra et al. 2012; Rit- 
ter et al. 2014). GRB progenitor models require 
massive, fast-rotating, and low-metallicity objects, 
(Galama et al. 1998; Woosley 1993; Hjorth et al. 
2003b; Woosley 2011; Berger et al. 2011; Levan 
et al. 2014), although the discovery of few GRBs 
that might have occurred in a solar or even su- 
per solar metallicity environment challenges this 
paradigm (e.g. Prochaska et al. 2009; Savaglio 
et al. 2012; Kriihler et al. 2012). 

Afterglow absorption spectroscopy of 2 ; > 1.5 
GRBs provides a unique tool to determine the con- 
stituents of the GRB environment, in particular 
the amount of metals produced by past and on- 
going star-formation in the vicinity of the GRB ex- 
plosion. This has important consequences for our 
understanding of the progenitor itself as well as the 
galaxies hosting GRBs throughout cosmic time. 
GRBs are identified, at first, based only on their 
high energy emission, therefore their hosts can be 
studied in great detail after the afterglow emission 
disappears, representing a sample of star-forming 
galaxies unbiased with respect to their intrinsic lu- 
minosity. In fact, it is possible to use GRBs and 
their hosts to trace cosmic star-formation in an 
independent way compared to magnitude limited 
Lyman-break galaxy surveys, although the effects 
of dust and metallicity biases are still under inves- 
tigation (Trenti et al. 2014; Bouwens et al. 2014; 
Levesque et al. 2010; Modjaz et al. 2008; Trenti 
et al. 2013; Kocevski & West 2011). 

Similar to GRBs, quasars (QSO) have also 
been used for decades to study the effects of re- 
ionization, the conversion of neutral hydrogen into 
stars, and the cosmic metal enrichment. In fact 
QSOs, like GRB optical afterglows, are very bright 
and can be seen up to very high redshift. The 
spectra of QSOs often show the presence of in- 
tervening absorbers (at redshifts lower than the 
QSOs), some of which are associated with large 
reservoirs of neutral hydrogen along their lines 
of sights. In particular. Damped Lyman-a sys- 
tems (DLAs), by definition, have column density 
of neutral hydrogen A^hi > 2xl0^°cm“^, while 
sub-DLAs are defined as absorbers with column 
density 10^® < A^m < 2xl0^°cm“^ (other types 
of subdivisions have been made, but they are not 
relevant for the purpose of this work). These ab- 
sorbers, which often trace galaxies along the QSO 


lines of sight (e.g. Fynbo et al. 2011; Schulze et al. 

2012) , are the best laboratories to investigate the 
ISM, its evolution, and cosmic star-formation at 
high redshift, providing important constraints on 
galaxy evolution models (e.g. Wolfe et al. 2005; 
Fynbo et al. 2011; Neeleman et al. 2013; Krogager 
et al. 2013; Jorgenson & Wolfe 2014; Ghristensen 
et al. 2014; Fumagalli et al. 2014; Rafelski et al. 
2012). 

Recently, Rafelski et al. (2014) have extended 
QSO-DLA studies up to z ~ 5: they show that 
the overall cosmological mean metallicity^ slowly 
decreases from 2 : « 1 up to z « 4.7 (see also Jor- 
genson et al. 2013; Prochaska et al. 2003; Rafelski 
et al. 2012) and then it appears to drop rapidly be- 
low the extrapolated linear metallicity evolution, 
as if a sudden metallicity enrichment in DLAs oc- 
curs shortly after the end of re-ionization. GRB- 
DLAs, in which the DLA is inside the GRB host, 
have been sparsely studied (Fynbo et al. 2008b; 
Savaglio et al. 2012; Sparre et al. 2013; Arab- 
salmani et al. 2014), mainly due to the small sam- 
ple size, the different data quality, and incomplete- 
ness. Nevertheless, Prochaska et al. (2007) have 
derived a higher metal content at z > 2 for a set 
of GRB-DLAs with respect to a large sample of 
QSO-DLAs, suggesting that GRBs probe denser, 
more dust depleted, and metal rich regions then 
the QSO-DLAs population (see also Fynbo et al. 

2013) . 

Finally, chemical evolution models suggest 
that DLAs metallicity measurements and rela- 
tive abundance ratios at very high-z enable us 
to better understand the effect of the primordial 
PopIII stars chemical enrichment (and IMF) onto 
subsequent PopII stars IMF (Salvador! & Ferrara 
2012; Kulkarni et al. 2013; Ritter et al. 2014). 

The goal of this study is twofold: 1) com- 

pare the findings by Rafelski et al. (2014), here- 
after R14, with a large sample of GRB-DLAs that 
extends previous studies (Prochaska et al. 2007; 
Savaglio 2012; Arabsalmani et al. 2014); 2) in- 
vestigate the metallicity evolution of GRB-DLAs 
and compare it with host galaxy metallicity pre- 
dictions at different redshifts. 

The paper is structured as follows: in the next 


^The cosmic mean metallicity is defined as (Z) =log 
(Xj lOlWifli JV(HI)i)/ Xj where i is the redshift bin 
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section we will present our samples of GRB- and 
QSO-DLAs; in §3 we will describe our analysis, in 
§4 we will discuss our findings and possible biases, 
and in §5 we will summarize our results. Through- 
out the paper we adopted the solar metallicity 
measurements from Asplund et al. (2009), with so- 
lar abundance of [X/H]q equal to 7.12 (Sulphur), 
7.51 (Silicon), and 4.56 (Zinc). 

2. Samples 

2.1. Gamma-ray Burst DLAs 

We select our GRB-DLAs sample from all the 
GRB afterglows observed during the 2000-2014 
time span for which HI and metallicity measure- 
ments can be obtained (Table 1). In order to de- 
tect the Lya absorption line (1216 A rest-frame) 
with most of the current spectrographs, a GRB 
has to be at least at z ^ 1.8 in order for the line 
to be redshifted out of the atmospheric blue cut-off 
(which usually means a minimum observed wave- 
length limit of ~3400 A). GRB-DLA absorbers 
are unambiguously associated with the GRB host 
galaxies, since often fine-structure transitions (e.g. 
Fe II* A2316) or the termination of Lya forest 
are identified at the same redshift of the Lya 
feature (Prochaska et al. 2006; Vreeswijk et al. 
2007). The presented sample includes spectra ob- 
tained with different instrument resolutions: from 
low resolution (resolving power R ^ 400) spec- 
tra obtained with the AIFOSC camera, to high 
resolution (i? ~ 55000) obtained with the UVES 
instrument. Due to the transient nature of GRB 
afterglows, some spectra were obtained when the 
afterglow was quite faint and the resulting signal- 
to-noise of the acquired spectra is not uniform 
within the sample. We therefore exclude GRBs 
with S/N < 3 (at usually 6000 A), due to the un- 
reliable detection of metal lines (see Section 3). 
All the data obtained with the Very Large Tele- 
scope (VLT) instruments were retrieved from the 
ESO Archive^ or were available within our collab- 
oration^. In few other cases we cannot obtain the 
raw data and we include the results from abun- 
dance analysis as they appear in the literature (e.g. 
GRB 050904 obtained with the FOCAS camera on 
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the Subaru telescope). We include these systems 
in Table I and briefly describe each line of sight in 
Appendix A. 

Our sample includes 13 GRB afterglow spectra 
obtained with the X-Shooter instrument mounted 
on the VLT. In order to analyze these data we 
primarily use our own customized pipeline writ- 
ten in IDL (Becker et ah, private comm.) as it is 
optimized for point sources and has an improved 
sky subtraction procedure. Additionally, we used 
the official pipeline (version 2.5, within the REFLEX 
workflow (Freudling et al. 2013; Goldoni et al. 
2006) to verify the output of our custom pipeline. 
For completeness, in the last column of Table 1 we 
report the literature reference where these data, if 
published, appear. 

Among the GRBs in Table 1 there are 12 
sub-DLAs and we exclude them from the sub- 
sequent analysis due to the fact that these ab- 
sorbers may probe a different environment than 
the general DLA population. These objects will 
be studied in a companion paper (Cucchiara et 
al. in prep), where we will present a more de- 
tailed description of the ionization held in order 
to reproduce the observed absorption pattern (see 
for example Vreeswijk et al. 2013). We also ex- 
clude 4 GRBs where only upper limits on the Vhi 
could be placed (GRB 071020, GRB 051111, and 
GRB 080913, GRB 090323), but we report them 
for completeness. For other 4 GRB-DLAs only 
Vhi measurement were obtained but the S/N is 
too low for reliably identifying any metal fea- 
ture (GRB 020124, GRB 060522, GRB080603B, 
and GRB121201A). In summary our GRB-DLA 
sample comprises 55 GRB-DLA lines of sight (3 
from the literature). Among these 55, we present 
metal abundances for 11 new GRB-DLA lines of 
sight. 

2.2. Quasar DLAs 

Thanks to the Sloan Digital Sky Survey large 
samples of quasars have been obtained up to 2 ^ 7. 
Thousands of these QSOs present DLAs along 
their lines of sight. We use the current most 
complete list of high-resolution QSO-DLA spectra 
obtained by Rafelski et al. (2012, hereafter R12) 
and R14, which extend previous work from, e. g., 
(Prochaska et al. 2003). All these QSOs have been 
observed with high-resolution spectrographs and 
have very high S/N, which are of great impor- 
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tance in order to resolve multiple narrow metal 
feature (see also Section 3.1) as well as to provide 
accurate metallicity measurements using different 
metal tracers. 

3. Analysis 

Our dataset allows us to assess ionic metal 
abundances directly from the afterglow spectra 
rather than relying on the literature. First, we 
select the spectra obtained with high-resolution 
spectrographs (which we defined at i? > 6000, 
see also Jorgenson et al. 2013) and high S/N, and 
compare them to the QSO-DLA sample that has 
been obtained with high-resolution instruments 
only (e.g. HIRES on the Keck telescopes). Sec- 
ond, as already presented by Prochaska (2006), 
abundance estimates (as well as the associate sta- 
tistical error) obtained using low-resolution spec- 
trographs are often under-estimated because of 
saturation effects, while blending may cause an 
over-estimation, and therefore the measurements 
should be cautiously used. As mentioned before, 
we perform our own analysis for the GRB line of 
sights and, in the following sections, we discuss in 
detail our results from the low-resolution sample. 
Our approach yields a more homogeneous sample, 
similar to the QSO-DLAs sample, for which the 
same procedure has been applied. 

3.1. Low Resolution Sample 

The consequences of low-resolution (R < 6000) 
spectroscopy on GRB-DLA abundance measure- 
ments via, e.g.. Curve of Growth methodology 
(COG, Spitzer 1978) has been already discussed 
in the literature (Prochaska 2006; Jorgenson et al. 
2013). Our large and diverse sample enables us 
to directly compare abundance measurements ob- 
tained for systems with different intrinsic metal- 
licity. In the first column of Figures 1 and 2 we 
show line profiles for two systems observed with 
the VLT /UVES instrument which provides a res- 
olution of 7 km s“^ (left column): GRB 050820 A 
(Figure 1) is a GRB-DLA with intrinsic moderate 
metallicity [ZjZ^ = —0.76), while GRB 050730 is 
a metal poor system {ZjZ^ = —1.96)). In the sec- 
ond and third column we present the same spec- 
tra resampled at the resolution of our average X- 
Shooter and Gemini/GMOS instruments respec- 
tively. Three important effects are evident and re- 


quire particular attention when estimating abun- 
dances with mid/low-resolution instruments: 1) 
blending of nearby lines, either from other tran- 
sitions rising from the GRB-DLA system or from 
unrelated intervening systems, may occur in lower 
resolution spectra (this also makes the determi- 
nation of the continuum level difficult, and even 
more so with low S/N data) 2) hidden saturated 
lines in low-resolution spectra may not be clearly 
identified and thereby yield a lower value then the 
true abundance; 3) strong absorption features as- 
sociated with moderately metal rich systems are 
still detected at i? ^ 1200 (or 200 km s“^, GMOS 
typical resolution) , but they completely disappear 
or are difficult to distinguish at even lower reso- 
lution, providing un-interesting metallicity abun- 
dance limits. 

In the first case, the result is an overestimation 
of metal abundances and de-blending procedures 
may be very complexed, especially when trying 
to asses further hidden saturation of the blended 
components. In the second, instead, the column 
density is underestimated. We therefore opted to 
use isolated weak lines (e.g. Fe 11I6O8), or very 
strong transitions (e.g. Si ii 1526) which are likely 
saturated and therefore provide reasonable lower 
limits. In the third case it is impossible to distin- 
guish between a true low-metallicity system and 
the effect of low-resolution instrumentation. 

Since the majority of our spectra (34 over 55) 
have been obtained with R < 2400 and the re- 
maining spectra with R > 6000 we adopt the lat- 
ter as the minimum resolution for which we can 
derive accurate metallicity estimates. Therefore, 
for all the spectra with resolution lower than the 
X-Shooter spectrograph (typically R < 6000) and 
S/N > 3, we measured metallicity from strong 
lines and provided only lower limits for the ionic 
abundances. We note that this analysis is simi- 
lar to Jorgenson et al. (2013), where ionic abun- 
dances from an even lower resolution instrument, 
like the MagE spectrograph on the Magellan tele- 
scope {R « 4200), was compared with the X- 
shooter, UVES and HIRES instruments. Also, 
while “hidden” saturation can be still present, we 
carefully choose features with depths such to min- 
imize this effect (Prochaska & Wolfe 1996; Pen- 
prase et al. 2010), typically with normalized flux 
values F\ < 0.5 in any pixel in the line profile. 
Jorgenson et al. (2013) have also preformed several 
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simulations on the reliability of using the Appar- 
ent Optical Depth method technique in deriving 
ionic abundances in such spectra (AOD, Savage 
& Sembach 1991) in comparison with Voigt pro- 
file fitting procedure. In order to further check 
this consistency we also performed a Voigt pro- 
file fitting for several ions in our X-Shooter and 
UVES spectra and compare these values with the 
ones obtained by our AOD analysis. The column 
densities agree with each other on average within 
one standard deviation. For these reasons (agree- 
ment between the AOD and Voigt profile meth- 
ods, the fact that the majority of our spectra have 
R < 6000, and that our comparison sample of 
QSO-DLAs metal abundances are also obtained 
with the AOD technique), we present our metallic- 
ity measurements derived with the AOD method- 
ology in Table 1. 

3.2. Neutral Hydrogen 

For all the GRBs in Table 1 for which we were 
able to retrieve the afterglow spectra, we deter- 
mine the redshift of the GRB-DLA based on the 
simultaneous identification of the strong Lya fea- 
ture (identifiable also at low resolution) and at 
least one of the fine-structure transitions (like 
Fell* and Nill*) often present in GRB afterglow 
spectra (Prochaska et al. 2006; Vreeswijk et al. 
2007; Prochaska et al. 2007). For those cases in 
which fine-structure lines were not observed, e.g. 
due to spectral coverage, we required, that at least 
other low ions transitions were detected at the 
same redshift of Lya, or that the end of the Lya 
forest was also identified. 

We fit the Lya profile with a Voigt profile using 
the x_fitdla procedure within the XIDL^ pack- 
age, while we adopted the measurements derived 
in the literature if the spectra were not available 
(see notes in Table 1). 

The QSO-DLAs measurements were obtained 
directly from R12 and R14, which performed a 
similar analysis on a sample of QSO spectra ob- 
tained with high-resolution spectrographs on the 
Keck telescope (e.g. HIRES or ESI). 

In Figure 3 we present the column density dis- 
tribution (fHi(N,z)) of our two samples in com- 
parison with the model by Noterdaeme et al. 
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(2009) (see also Wolfe et al. 1995; Peroux et al. 
2003; Prochaska & Wolfe 2009; Noterdaeme et al. 
2012b). As previously noted (Reichart & Price 
2002; Savaglio et al. 2003; Jakobsson et al. 2006; 
Prochaska et al. 2007), HI column density in QSO- 
DLAs is a factor of ten lower than the GRB-DLAs 
(although see Noterdaeme et al. 2012a, 2014) . 
This suggests that GRB-DLAs may trace a denser 
ISM phase, more similar to the sites of ongoing 
star formation (Fynbo et al. 2008b), while QSO- 
DLAs may probe, instead, a lower column density 
medium, possibly further from dense molecular re- 
gions. Moreover, a handful of GRBs with large 
amount of neutral hydrogen (Ahi > 10^^'®cm“^) 
along their line of sights exhibit the presence of 
molecular hydrogen (H2) (Prochaska et al. 2009; 
Kriihler et al. 2013; D’Elia et al. 2014), while only 
in a few cases have been found along quasars (No- 
terdaeme et al. 2008; Srianand et al. 2012; Jorgen- 
son et al. 2013, 2014). 

3.3. Metallicity 

Measuring the gas metal content in GRB-DLA 
systems is not simple, especially because some of 
the metals might be locked into dust grains (dust 
depletion effect, Savaglio et al. 2003). For exam- 
ple, the mildly refractory element silicon is usu- 
ally depleted in the Galaxy, but only marginally in 
DLAs and therefore it can often be used to deter- 
mine the gas metallicity (Wolfe et al. 2005; Rafel- 
ski et al. 2012), assuming that this also applies for 
GRB-DLAs. Other good metallicity tracers in- 
clude sulfur and zinc. Zinc, in particular, is often 
preferred because it is un-depleted in the ISM and 
has two strong transitions at rest-frame 2026 A 
and 2063 A. However, it is only a trace element and 
therefore represents a small fraction of the mass 
density of the heavy element. Moreover, the evo- 
lution of zinc resembles iron only for specific star 
formation histories and careful modeling of zinc 
production by SN II and SN la shows an under- 
production of Zn and Mg compared to S, invoking 
additional production site, such as intermediate- 
mass stars, in order to reconcile their abundance 
values(see Fenner et al. 2004, for a detailed descrit- 
pion). Therefore, as pointed out by Rafelski et al. 
(2012) and Prochaska (2006) we decided to use 
low-ionized transitions of sulfur (e.g. S Iil250), 
silicon (e.g. Si 11I8O8), and iron (e.g. Fe 11I6II), 
in order of importance. We only use the zinc lines 
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in the few cases where the previously discussed 
lines are unavailable, such as when they fall at the 
location of atmospheric telluric bands. 

We measured ionic abundances using the Ap- 
parent Optical Depth method (Savage & Sembach 
1991), which relies on the identification of several 
unsaturated lines of the same species and provides 
accurate measurements with no assumptions on 
the features Doppler parameter as other methods 
(like Curve-of-Growth analyses, Carroll & Ostlie 
1996). As we mention in the previous earlier these 
values well agree with the Voigt profile fitting tech- 
nique. 

In some cases, due to the presence of saturated 
transition, we only estimate lower limits on the 
overall metallicity (see also Savaglio & Fall 2004; 
Savaglio 2006; Prochaska 2006). In fact, to be con- 
servative we consider all our measurements from 
low-resolution data as lower limits. 

In Table 1 we summarize our findings, includ- 
ing HI and metallicity measurements, as \X/H] 
relative to solar^, and the ion used along that spe- 
cihc line of sight. No dust depletion or ionization 
correction have been applied (see Sections 4.3 and 
4.2). In the last column we list the references rel- 
ative to each GRB and, in the Appendix A, we 
briefly described each line of sight in more detail. 

4. Discussion 

4.1. Possible Observational Biases 

First of all we need to understand if there are 
possible biases that might affect either samples: 
GRBs are selected solely on the prompt emission 
detection in the gamma-ray energy bands from 
dedicated satellites (e.g. Swift). The spectra were 
taken mostly independently of the brightness of 
the afterglow, although a preference in observing 
brighter events might be present. Neither of these 
two selection criteria (brightness in the gamma-ray 
bands or in the optical) seem to explain the differ- 
ence in TVhi distribution (Figure 3). As previously 
noted by Pontzen et al. (2010), GRB-DLAs seem 
to probe a different type of absorbers than the 
QSO-DLA population. In particular, if the faint 
afterglows were not observed spectroscopically be- 
cause some high level of extinction (like the “dark” 
GRBs from Perley et al. 2013), and assuming a 


^[X/H] = log[X/H]oLA-log[X/H]Q 


possible correlation between metallicity and visual 
extinction (see Zafar & Watson 2013), this would 
bias us towards lower metallicities. This, com- 
bined with a small number of metal-poor GRB- 
DLAs at high- 2 ; suggests that, if such bias exists 
it would have a small effect, also because at high 
redshift there are less dusty systems (Govino et al. 
2013). At 2 > 6, our understanding of dust pro- 
duction mechanisms (e.g. SN, AGB stars) are still 
to be fully understood, although some advance- 
ments have been made both theoretically and ob- 
servationally (Gall et al. 2011). We also note that, 
similarly to these authors, we are assuming the ex- 
tinction along the line of sight as the same as the 
one derived from afterglow studies, which in prin- 
ciple can lead to different estimates with respect 
to the hosts galaxy extinction (see Elliott et al. 
2013; Perley et al. 2013, for a comparison between 
the afterglow- and host-derived extinction) . Over- 
all, despite the large sample presented here, our 
conclusions, in particular at high redshift, require 
more data in order to better assess the effect of 
such observational biases. 

The QSOs sample from R12 and R14, instead, 
has been selected solely on the presence of a Lya 
line in their SDSS spectrum and logVni > 20.3, 
therefore represent an unbiased sample with re- 
spect to the gas metallicity (Rafelski et al. 2012). 
The z <2 QSO-DLAs were selected based on the 
presence of Mg ii in the SDSS spectra, implying a 
small bias against low-metallicity systems at such 
redshift. Since our redshift range of interest is 
mostly at 2 : > 2, this should not be a source of 
significant concern. 

4.2. Ionization Correction 

The identification of fine-structure transitions 
and the line profile variability within the first 
hour after the GRB explosion have been associ- 
ated to the effect of an evolving radiation field 
from the GRB (like the afterglow) or nearby young 
stars onto the progenitor surrounding medium 
(Prochaska et al. 2006; Vreeswijk et al. 2007; De 
Gia et al. 2012). This processes directly affects 
metal abundance measurements and the inferred 
metallicity ionization correction is strongly depen- 
dent on the hydrogen column density (Vreeswijk 
et al. 2013). 

Distinguishing between ionized (photoionized 
or excited) gas from the GRB emission in the 
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vicinity of the GRB progenitor and the bulk of the 
host galaxy ISM, which is also ionized by the sur- 
roundings radiation field, is a complicated process: 
if the gas probed by the afterglow spectroscopy is 
close to the GRB (< 200pc), then the ionization is 
increasing with time and complex modeling is re- 
quired to estimate the time dependent ionization 
correction (Krongold & Prochaska 2013). Proba- 
bly, the most important example in which a de- 
tailed photoionization modeling (including photo- 
ionization and excitation effects) has been per- 
formed, is the sub-DLAs GRB 080310 (Vreeswijk 
et al. 2013). In most of our GRB-DLAs the data 
have been acquired at much later time and in sev- 
eral cases, when multiple spectra were obtained, 
these objects do not present line profile variation: 
such non-variability suggests that we are likely 
probing the ISM at larger distances and/or gas 
clouds unaffected by the GRB radiation. 

Reassuringly, similar models to the ones of 
Vreeswijk et al. (2013) have been performed in 
GRB-DLAs spectra (Ledoux et al. 2009): ion- 
ization corrections are often minimal (< 10%) 
compared to other even extreme cases, e.g. star- 
forming galaxies or Lyman-a emitters. For these 
reasons we did not apply any ionization correc- 
tion to our derived column densities and our final 
metallicity measurements. We stress that the ion- 
ization correction is an important aspect of the 
column density determination and requires par- 
ticular attention and a much larger sample of 
rapid high-resolution spectroscopic observation se- 
quences starting not later than a few minutes after 
the burst. While this is the least known quantity 
in our study we argue that our results are indica- 
tive of an overall general characterization of the 
GRB-DLA population. A large sample of multi- 
epochs high-resolution datasets and accurate mod- 
eling is still needed, but is beyond the scope of this 
paper. 

4.3. Depletion Correction 

Understanding the effect of dust depletion in 
DLAs requires the determination of ionic abun- 
dances of both refractory and non-refractory el- 
ements (see for example D’Elia et al. 2014). 
Savaglio & Fall (2004) determined for the first time 
such “depletion pattern”, while recently De Cia 
et al. (2013) compared a sample of 20 GRB-DLAs 
with 47 QSO-DLAs in order to study the dust-to- 


metal ratio of these systems. These authors used 
the [Fe/Zn] as dust indicator and found that the 
dust depletion correction is in the most depleted 
cases ^ -1-0.1 dex independently of metallicity. 

These results are in contrast with R12 where 
strong iron depletion is present at ZJZq > —1.0. 
Instead, R12 shows that silicon is rarely depleted 
in QSO-DLAs and that depletion effects are rel- 
evant only at high metallicity (> — O.S.Z/Zq), 
therefore are minimal for the majority of our GRB 
lines of sights(see also Vladilo et al. 2011). 

For our GRB-DLA sample, to calculate a deple- 
tion correction is not always possible, especially 
for our low resolution GRB-DLA spectra. Also, 
in the few cases where iron was used, we derived 
ZjZQ < —1.0, so depletion correction should be 
minimal. 

Nevertheless, the application of a depletion cor- 
rection would increase the metallicity of our sys- 
tems, strengthening some of our conclusions (see 
Sections). Therefore, we choose not to apply 
any depletion correction to the measured metal- 
licities. Finally, we performed a-enhancement 
correction to the iron based metallicities for our 
high-resolution sample using the correction factors 
adopted by R12. 

4.4. Metallicity Evolution 

Our hndings are shown in Figure 4, where the 
QSO-DLAs (in gray) and GRB-DLAs (in red) 
metallicities are plotted. We perform a linear 
fit of the metallicity measurements (and lim- 
its) with redshift using a survival analysis tech- 
nique (Schmitt 1985), which takes into account 
upper and lower limits. In particular we used 
the statistics . schmittbin package within the 
IRAF® distribution. We also used a bootstrap sam- 
pling in order to determine the 1-sigma error in 
the fit (with 500 iterations). 

Ideally, we would like to use our metallicity 
measurements to investigate the cosmic metal 
budget at different epochs: this is usually done 
weighting the average metallicity over a specific 
redshift bin with the total neutral hydrogen col- 
umn density in the same redshift interval (see 


®IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation. 
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R12 and R14 for the QSO-DLAs sample). Un- 
fortunately, because the large number of limits 
present in the GRB-DLA sample we simply fit the 
metallicity of the single systems, which still pro- 
vide useful insight on the DLA populations metal 
content. 

For the GRB-DLA sample, we derive [X/ H]grb = 
(—0.07 ± 0.06)z — (0.75 ± 0.25) (thick red dashed 
line), which is also consistent with no-evolution 
(at the Icr confidence level). For the QSO-DLA 
sample we derive a linear trend between red- 
shift and DLA metallicity given by [X/H] = 

(— 0.20 ±0.03)z — (0.68 ±0.09) (thick black dashed 
line in Figure 4). 

To test the reliability of our fit and how much 
this is sensitive to the small number of abundance 
measurements (with respecting the limits) we ran 
a series of Monte Carlo simulations: we created 
1000 mock samples of GRB-DLAs metallicity mea- 
surements and limits (with values within the typ- 
ical GRB-DLA high-resolution points RMS) and 
we repeated the survival analysis fit. Assuming an 
intrinsic true slope from the QSO-DLA distribu- 
tion of —0.2 we obtained a slope > —0.07 or less 
only in 7% of the cases, which means that it is 
unlikely that our results are affected by low num- 
ber statistics. Similar tests with different intrin- 
sic distribution (from flat to very steep metallicity 
evolution) provide similar results and reassure us 
that, despite the small number statistics, we can 
recover the input slopes to within the reported Icr 
confidence interval. 

Finally, we point out that GRB-DLAs and 
QSO-DLAs have a similar metallicity distri- 
bution at z ^ 2 — 3, suggesting that there 
is no difference among the two populations 
of absorbers in terms of metal content. We 
performed a two-sample Kolmogorov-Smirnov 
test among these subsamples using the IRAF 
stsdas . analysis . statistics .twosampt task 
and we can rule out the null hypothesis that the 
two distributions are drawn from the same parent 
population at the 90% confidence level. Instead, 
if we naively consider the lower limits as actual 
measurements and apply an arbitrary correction of 
different values (from ±0.1 to ±0.3 dex), there re- 
mains evidence that the two distributions are dis- 
tinct even with the highest value considered (±0.3 
dex), but the statistical significance is < 99%. 

However, the GRB-DLA metallicity declines 


suggests that, in particular at z > 3, GRB-DLA 
environment is more metal enriched than in QSO- 
DLAs, likely by active star-formation episodes. 
These metals may have been ejected by super- 
nova explosions or mass losses and polluted the 
GRB progenitor’s neighborhood before the GRB 
occurred (Kroger et al. 2006; Mao 2010; Matteucci 
& Recchi 2001; Kulkarni et al. 2013). If this is 
true, the inferred metallicity may not reflect the 
overall metal content of such high -2 GRB host 
galaxies, although we know that the gas inter- 
cepted by GRB afterglow spectra may lie up to 
few kpc from the GRB explosion site (see next 
section). 

4.5. Characterizing these DLA popula- 
tions 

We will now try to understand the different 
metallicity evolution between the GRB- and QSO- 
DLA samples. For example, if GRB-DLAs trace 
the general host galaxies ISM, this gas ionization 
state may show the effect of whatever local ion- 
ization field is in the vicinity of the GRB or along 
the line of sight. 

Prochaska et al. (2007) argued that the size 
of such molecular cloud would exceed the largest 
molecular cloud in the local group, so the GRB- 
DLAs are tracing material as far as at least lOOpc 
from the GRB (and even out to ~ 2 kpc in the case 
of GRB 060418, Vreeswijk et al. 2007). The fact 
that in the z = 3.5 — 6 redshift range the metallic- 
ity of the gas is, on average 10% of the solar value 
suggests that a substantial amount of metals are 
already present at high redshift. 

These estimates are in agreement with the most 
recent GRB progenitor models (Woosley 2011; 
Woosley & Heger 2012): if the GRB host galaxy 
is metal rich, then a large amount of metals have 
been produced by SNe or, in the highest redshift 
bin observed, by a late population of PopIII stars 
(or more likely early PopII). Nevertheless, in or- 
der to be able to produce a GRB explosion from 
Wolf-Rayet type stars, the amount of metal in- 
jected into the ISM and mixed throughout the 
host galaxy has to be below a certain threshold to 
retain enough angular momentum and minimize 
the stellar mass loss of the progenitor (Woosley 
& Heger 2006; Woosley 2011; Woosley & Heger 
2012). Clearly more data need to be acquired, 
in particular at high -2 to confirm the existence 


of such limits, also in light of the high-metallicity 
hosts observed (e.g. Elliott et al. 2013; Kriihler 
et al. 2012; Savaglio et al. 2012). 

On the other hand, DLAs identified in QSO 
spectra are cross-section dependent, and there- 
fore they might not trace the denser part of the 
galaxies (like for GRB-DLAs). Also, it has been 
proposed that a combined large sample of GRB- 
and QSO-DLAs may represent a complete cen- 
sus of z Ri 3 star-forming galaxies that could be 
missed by magnitude limited surveys (Fynbo et al. 
2008b). 

Finally, as R14 pointed out, the observed metal- 
licity decrement at z > 4.7 suggests an increase 
in the covering fraction of neutral gas: similar 
behavior in fact can be produced by the combi- 
nation of increased density of the Universe and 
lower background radiation field, which allows the 
hydrogen gas to be self-shielded at the lower den- 
sity as function of redshift (Fumagalli et al. 2013). 
In this picture, the denser region would reside in 
the halo of the galaxies or in the IGM where the 
star- format ion, and therefore metal enrichment, is 
lower. 

4.6. GRB-DLA metallicity in context 

A long standing debate exists as to the degree to 
which GRBs faithfully trace the cosmic star for- 
mation rate. Although GRBs have been associ- 
ated with broad lined supernovae at low redshift 
and regions of active star formation in their host 
galaxies, spectroscopic observations have shown 
that GRB host galaxies tend to be relatively metal 
poor compared to SNe Ibc hosts (Modjaz et al. 
2008). This observed preference for low metallic- 
ity environments may impart a redshift dependent 
bias in the type of star forming regions that can 
produce a GRB (Kocevski et al. 2009; Trent! et al. 
2014). 

In particular, at low redshifts (z < 1), a prefer- 
ence for low metallicity environments would limit 
GRBs to low mass spirals and dwarf galaxies (e.g. 
Levesque et al. 2010), due to the well established 
relationship between mass and metallicity (see Ko- 
cevski & West 2011; Graham & Fruchter 2013, 
and references therein). At higher redshifts, the 
mass range of galaxies capable of hosting a GRB 
would increase to include more massive, star form- 
ing galaxies, since the average metallicity of all 


galaxies in the Universe falls. Recent unbiased 
searches of GRB host galaxies like the THOUGH 
survey (Hjorth et al. 2012) or of the host galaxies 
of the “dark” GRB population (Perley et al. 2013) 
largely support this trend. These surveys find that 
bursts at intermediate redshifts tend to be drawn 
from star forming galaxies with a greater diversity 
of mass, morphology, and dust content, suggest- 
ing that high redshift GRBs may serve as more 
faithful tracers of cosmic star- format ion compared 
to their low redshift counterparts (see also Hunt 
et al. 2014; Perley et al. 2014). 

Our sample of GRB-DLAs, although not a com- 
plete host sample, covers a much greater redshift 
range than the emission line derived measurements 
in these studies and, more importantly, does not 
depend on strong observational biases (e.g. bright- 
ness of the host, intensity of the emission lines), 
although we may miss some dusty (therefore metal 
rich) events. Based on the arguments outlined 
above, even as metallicity biased tracers of star 
formation, the GRB-DLA results presented in Fig- 
ure 4 should become more representative of the 
metallicity evolution of the general star forming 
galaxy population with increasing redshift. How- 
ever, a better understanding of both observational 
biases (Fynbo et al. 2008a, b) and the effect of 
metallicity in GRB production (Jakobsson et al. 
2013) are required before we can fully address the 
connection between the environments which are 
capable of producing GRBs and the conditions of 
star forming regions in the early Universe. 

From a theoretical standpoint, recent simula- 
tions by Trent! et al. (2014) have been suggest- 
ing that two combined channels of GRB popu- 
lations exist: one, the “collapsar” mode, which 
strictly depends on the host metallicity, and a sec- 
ond one, the “binary stars” mode, which is metal- 
licity independent. Assuming these two modes 
coexist and the known mass-metallicity relation 
of star forming galaxies, these authors predict 
GRB host galaxies metallicity redshift evolution 
with different combinations of these two channels 
(from strong metallicity bias to an almost negligi- 
ble one). Unfortunately, metallicity measurements 
from emission line diagnostics are not yet available 
for z > 1 host galaxies, but our sample represents 
the best opportunity to test these models. 

In Figure 5 we present our metallicity results in 
comparison the predicted metallicity for the up- 
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per 95%, the median, and the bottom 20% of the 
GRB host galaxies distribution, assuming an ab- 
sent (dotted lines) and a moderate (solid lines) 
metallicity bias (adapted from Trenti et al. 2014, 
using a value of p = 0.04). 

These models essentially predict that in any 
given redshift bin, for example, 20% of the hosts 
that produce a GRB have metallicity below the 
lower solid line. Indeed, focusing on the redshift 
z = 2 — 3 range, 5 over 27 hosts have metallic- 
ity below the line. On the other hand the “no- 
metallicity bias” model does not agree with the 
data. Overall, it seems that at least at z ;< 4 a 
moderate metallicity bias is required in order to 
reconcile theory and observations (Vergani et al. 
2014), although a more detailed analysis of these 
models and the implications on the metallicity cut- 
off in the GRB host masses is needed, in particular 
to understand the role of such metallicity bias, if 
indeed exists, at higher redshift. 

5. Conclusions 

In this work we investigated the properties of 
GRB lines of sight that show evidence for Damped 
Lyman-a systems within the GRBs host galaxy. 
Similar systems have been studied along QSOs in 
order to understand the chemical enrichment and 
the metal evolution in cold gas systems over cos- 
mic times. 

We collected all the publicly available GRB af- 
terglow spectra, including data already published 
in the literature and we uniformly analyzed these 
lines of sight in order to obtain metal abundances 
using mildly or non refractory species (e.g. sul- 
phur, silicon, and zinc). In particular, we opted 
for a conservative approach when considering low- 
resolution data (see Section 3.1): we consider de- 
pletion correction to our metal abundance mea- 
surements a negligible effect (see Section 4.3). 
Also, we note that ionization correction may have 
an important role in estimates derived from very 
early time spectroscopic data (< 10 minutes), due 
to the highly variable ionization flux from the 
GRB itself. Nevertheless, since most of our data 
have been acquired at much later times and the 
fact that detailed analyses are not always possi- 
ble (especially with low-resolution data) we do not 
consider such correction in our analysis (see Sec- 
tion 4.2). 


In Figure 4 we present our GRB-DLAs metal- 
licity evolution ht in comparison with the QSO- 
DLA results from R12. We performed a detailed 
survival analysis fit to accurately take into ac- 
count all the limits present in our dataset. As in 
previous works (Savaglio et al. 2012; Arabsalmani 
et al. 2014) we derive a much shallower decline 
of the GRB host metallicity with redshift rela- 
tive to the QSO-DLA sample suggesting a some- 
what metal rich environment for the GRB host 
galaxies (in particular at z > 4), but still be- 
low solar values, in agreement with GRB progen- 
itor models (Woosley & Heger 2012). At z < 3 
there is a reasonable overlap between the two 
populations, according to their metallicity proper- 
ties, indicating similarities between the two DLA 
populations, despite the higher HI column den- 
sity traced by the GRB-DLAs. At higher red- 
shift, GRB-DLAs seem to prefer higher metallic- 
ities than QSO-DLAs. Despite the small number 
of high-resolution data, a few lines of sight (e.g. 
GRB130606A or GRB060510B) seem to indicate 
that these GRB-DLAs point towards a denser, 
metal rich environment, likely tracing a less com- 
mon population of metal-rich DLAs. 

Overall, our findings confirm the idea that mov- 
ing towards higher redshifts, GRBs trace prefer- 
entially a denser metal rich environment within 
galaxies, while the QSO-DLA population may be 
progressively dominated by neutral regions with 
minor star formation. In other words, the ab- 
sence of metal poor {[X/ H]crb ^ —1.5) GRB- 
DLAs at high-z seems to indicate that such GRB- 
DLA hosts are at the low-end of the luminosity 
(or mass) distribution and do not present high 
star- format ion (and GRB progenitors). Another 
possibility is that these GRBs occurred in a com- 
pletely different environments (galactic haloes). 
Deep imaging, with 8-10 meter telescopes or with 
HST, is needed in order to fully characterize these 
hosts. 

Finally, we compare our findings with the 
most recent model prediction of GRB host galaxy 
metallidties (Figure 5, Trenti et al. 2014). Our 
results broadly agree with host galaxy metallicity 
predictions where two-channels for GRB produc- 
tions are considered (collapsars and binary) and 
mildly affected by a metallicity bias. The model 
predicts well the metallicity distribution for the 
bottom 20% and median of the host population. 
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though more work is required to fully understand 
this metallicity bias and its effect on other prop- 
erties of the GRB host population (e.g. the hosts 
mass limit driven by a metallicity cut-off Kocevski 
et al. 2009). 

Our GRB-DLA host galaxies represent the 
largest sample available to date and, although not 
complete, it is suitable for multi band follow-up, 
in particular for the current and upcoming near- 
infrared spectroscopic instruments, which will al- 
low to determine SFR and metallicities directly 
using emission line diagnostics. Also, the deter- 
mination of the hosts properties, like mass and 
rest-frame UV star formation, will help in bet- 
ter characterizing the overall high-z GRB host 
population, and their capability to harbor GRB 
progenitors. Furthermore, we will be able to bet- 
ter understand the observational biases that might 
affect our results, especially at high-z, where few 
lines of sights are observed (including dust extinc- 
tion and/or cosmic metallicity trend). 

While multiband surveys of this sample of 
GRB-DLA hosts will allow a better characteri- 
zation of these galaxies, the advent of future mis- 
sions like JWST and the new generation of 30-m 
telescopes will be able to identify these faint hosts 
at the highest redshifts and spatially resolve the 
regions of star-formation traced by GRB-DLAs, 
which seem to hide the secrets of primordial star 
formation sites. 
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Table 1 

GRB-DLAs Sample 


GRB 

ZGRB 

log(AfHi) 

[X/H]- 

Ion 

Fine- 

Structure 

Telescope/ 

Instrument 

Resolution 

at 6000A 

S/N 

per pixel 

Reference 

000926 

2.3621 

21.3 ±0.25 

> -0.30 

Zn 

N 

Keck/ESI 

20000 

10 

[1] 

011211 

2.1427 

20.4 ±0.2 

> -1.22 

Si 

N 

VLT/FORS2 

2400 

10 

[2] 

020124 

3.198 

21.7 ±0.2 



N 

VLT/FORSl 

450 

4 

[3] 

021004 

2.3289 

19.0 ±0.2 



Y 

VLT/UVES 

40000 

6 

[4] [5] 

030226 

1.98 

20.5 ± 0.3 

> -1.28 

Fe 

Y 

Keck/ESI 

20000 

40 

[6] 

030323 

3.3714 

21.9 ±0.07 

> -1.32 

S 

Y 

VLT/FORS2 

2100 

20 

[7] 

030429 

2.658 

21.6 ±0.2 

> -1.13 

Si 

N 

VLT/FORSl 

600 

40 

[8] 

050319 

3.24 

20.9 ±0.2 

> -0.77'= 

S 

N 

NOT/AIFOSC 

355 

4 

[9] [10] 

050401 

2.899 

22.6 ±0.3 

> -1.07 

Zn 

Y 

VLT/FORS2 

545 

10 

[11] 

050505 

4.27 

22.05 ±0.1 

> -1-2 

S 

Y 

Keck/LRIS 

1200 

20 

[12] 

050730 

3.96723 

22.1 ±0.1 

-1.96 ±0.11 

s 

Y 

VLT/UVES 

40000 

10 

[13] [14] 

050820A 

2.6145 

21.1 ± 0.1 

-0.76 ±0.13 

s 

Y 

VLT/UVES 

40000 

12 

[8] 




-0.78 ±0.11 

Fe 

Y 

Keck/HIRES 

30000 

10 

[14] 

050904 

6.26 

21.3 ±0.2 

> -1.0'= 

S 

Y 

Subam/FOCAS 

1000 

7 

[15] 

050908 

3.344 

19.4 ±0.2 



N 

Gemini/GMOS 

1200 

20 

[16] 

050922C 

2.1996 

21.55 ±0.1 

-1.88 ±0.14 

s 

Y 

VLT/UVES 

45000 

10 

[17] 

051111 

1.549 

< 21.9 



Y 

Keck/HIRES 

55000 

20 

[18] 

060115 

3.533 

21.5 ±0.1 

> -1.53 

s 

Y 

VLT/FORSl 

990 

4 

[9] 

060124 

2.3 

18.5 ± 0.5 



N 

Keck/LRIS 

1200 

18 

[9] 

060206 

4.048 

20.85 ±0.1 

> -0.74 

s 

Y 

Lick/KAST 

1200 

28 

[9] 

060210 

3.913 

21.55 ±0.15 

> -0.83 

Si 

Y 

Gemini/GMOS 

1200 

40 

[9] 

060223A 

4.41 

21.6 ±0.1 

> -1.8'= 

s 

N 

Keck/LRIS 

1200 


[19] 

060510B 

4.94 

21.3 ±0.1 

> -0.84 

s 

N 

Gemini/GMOS 

1200 

15 

[19] 

060522 

5.11 

21.0 ±0.3'= 



N 

Keck/LRIS 

1200 

2 

[19] 

060526 

3.221 

19.9 ±0.15 



N 

VLT/FORSl 

1200*’ 

18 

[20] 

060605 

3.773 

18.9 ±0.4 



Y 

PMAS 

500 

7 

[21] 

060607A 

3.075 

16.95 ±0.03 



Y 

VLT/UVES 

55000 

30 

[17] 

060707 

3.425 

21.0 ±0.2 

> -1.69 

Fe 

Y 

VLT/FORS2 

800 

7 

[22] 

060714 

2.711 

21.8 ±0.1 

> -0.97 

Zn 

Y 

VLT/FORSl 

800 

30 

[22] 

060906 

3.686 

21.85 ±0.10 

> -1.72 

S 

N 

VLT/FORSl 

800 

8 

[22] 

060926 

3.206 

22.6 ±0.15 

> -1.32 

Zn 

Y 

VLT/FORSl 

800 

20 

[22] 

060927 

5.464 

22.50 ±0.15 

> -1.55 

S 

N 

VLT/FORSl 

500 

3 

[23] 

061110B 

3.433 

22.35 ±0.10 

> -1.84 

S 

Y 

VLT/FORSl 

800 

10 

[9] 

070110 

2.351 

21.7 ±0.1 

> -1.32 

Si 

Y 

VLT/FORS2 

800 

15 

[9] 

070411 

2.954 

19.3 ±0.3 



Y 

VLT/FORS2 

800 

8 

[9] 

070506 

2.308 

22.0 ±0.3 

> -0.65 

Zn 

N 

VLT/FORSl 

800 

18 

[9] 

070721B 

3.628 

21.5 ±0.2 

> -2.14 

Si 

Y 

VLT/FORS2 

800 

5 

[9] 

070802 

2.455 

21.5 ±0.2 

> -0.54 

Si 

Y 

VLT/FORS2 

800 

7 

[9] 

071020 

2.145 

< 20.30 



N 

VLT/FORS2 

800 

5 

[9] 

071031 

2.692 

22.15 ±0.05 

-1.85 ±0.12 

Fe 

Y 

VLT/UVES 

55000 

10 

[24] 

080210 

2.641 

21.9 ±0.1 

> -1.37 

Fe 

Y 

VLT/FORS2 

1400 

25 

[9] [25] 

080310 

2.427 

18.7 ±0.1 



Y 

VLT/UVES 

55000 

30 

[17] 

080413A 

2.433 

21.85 ±0.15 

> -1.56 

Zn 

N 

Gemini/GMOS 

1200 

17 

[24] 

080603B 

2.69 

21.85 ±0.05^= 



Y 

NOT/AIFOSC 

355 


[9] 

080607 

3.037 

22.7 ± 0.15 

> -1.72 

Fe 

Y 

Keck/LRIS 

2000 

40 

[9] [26] 

080721 

2.591 

21.6 ±0.1 

> -1.73 

S 

N 

VLT/FORSl 

800 

40 

[27] 

080804 

2.20542 

21.3 ±0.1 

-0.75 ±0.16 

Zn 

N 

VLT/UVES 

55000 

10 

[9] 




^ -1-25 

Zn 

N 

Gemini 

1200 

10 

[10] 

080810 

3.35 

17.5 ±0.15 



Y 

NOT/AIFOSC 

400 

8 

[28] 






Y 

Keck/HIRES 

50000 

30 

[28] 
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Table 1 — Continued 


GRB 

ZGRB 

log(AfHi) 

[X/H]- 

Ion 

Fine- 

Structure 

Telescope/ 

Instrument 

Resolution 

at 6000A 

S/N 

per pixel 

Reference 

080913 

6.69 

< 19.84 



N 

VLT/FORS2 

800 

3 

[29] 

081008 

1.96 

21.59 ±0.1 

-0.86 ±0.14 

s 

Y 

VLT/UVES 

40000 

5 

[30] 




> -1.41 

s 

Y 

Gemini/GMOS 

1200 

12 

this work 

090205 

4.64 

20.73 ±0.05 

> -0.57 

s 

Y 

VLT/FORSl 

440 

5 

[31] 

090323 

3.5778 

> 19.90 



Y 

VLT/FORS2 

1200 

24 

[32] [33] 

090426 

2.609 

19.1 ± 0.15 



N 

Keck/LRIS 

1200 

7 

[20] 

090516 

4.109 

21.73 ±0.1 

> -1.36 

Si 

Y 

VLT/FORS2 

800 

70 

[34] 

090809 

2.73 

21.40 ±0.08 

-0.57 ±0.10 

Si 

Y 

VLT/X-Shooter 

8000 

12 

[35] 

090812 

2.425 

22.3 ±0.1 

> -1.64 

Si 

Y 

VLT/FORS2 

800 

60 

[36] 

090926A 

2.1062 

21.73 ±0.07 

Ri -1.9 

s 

Y 

VLT/FORS2 

780 

35 

[37] 




-2.18 ±0.12 

s 

Y 

VLT/X-Shooter 

10000 

20 

[37] 

100219A 

4.667 

21.13 ±0.12 

-0.95 ±0.18 

s 

Y 

VLT/X-Shooter 

6000 

6 

[38] 




IV 

1 

bo 

s 

Y 

Gemini/GMOS 

1200 

7 

this work 

100425A 

1.755 

21.05 ±0.10 

-0.96 ±0.42 

Fe 

Y 

VLT/X-Shooter 

8000 

4 

[35] 

110205A 

2.214 

21.45 ±0.2 

> -0.82 

s 

Y 

FAST 

2400 

30 

[39] 

111008A 

4.98968 

22.3 ±0.06 

-1.63 ±0.13 

Fe 

Y 

VLT/X-Shooter 

10000 

10 

[40] 

111107A 

2.893 

21.0 ±0.10 

> -0.45 

S 

Y 

VLT/X-Shooter 

8000 

5 

this work 

120327A 

2.813 

22.01 ±0.09 

-1.51 ±0.11 

S 

Y 

VLT/X-Shooter 

8000 

30 

[41] 

120716A 

2.487 

21.55 ±0.15 

> -1.76 

Fe 

Y 

VLT/X-Shooter 

8000 

7 

this work 

120815A 

2.3574 

21.95 ±0.1 

-0.93 ±0.13 

Zn 

Y 

VLT/XShooter 

10000 

12 

[42] 

120909A 

3.9293 

21.20 ±0.10 

-0.66 ±0.11 

S 

Y 

VLT/X-Shooter 

8000 

9 

this work 

121024A 

2.2977 

21.50 ±0.10 

-0.40 ±0.12 

Zn 

Y 

VLT/X-Shooter 

8000 

15 

[43] 

121201A 

3.385 

21.7 ± 0.2 




VLT/X-Shooter 

8000 

;S 3 

this work 

130408A 

3.757 

21.70 ±0.10 

-1.24 ±0.12 

S 

Y 

VLT/XShooter 

8000 

50 

this work 




> -1.1 

s 

Y 

Gemini/GMOS 

1200 

20 

this work 

130505A 

2.2687 

20.65 ±0.10 

> -1.42 

Fe 

Y 

Gemini/GMOS 

1200 

30 

this work 

130606A 

5.9134 

19.93 ±0.2 



Y 

Gemini/GMOS 

1200 


[44] 

140226A 

1.9733 

20.60 ±0.20 

> -0.54 

Fe 

N 

keck/LRIS 

1200 

30 

this work 

140311A 

4.953 

21.80 ±0.30 

> -1.65 

Ni 

Y 

Gemini/GMOS 

1200 

7 

this work 

140419A 

3.961 

19.3 ±0.2 



Y 

Gemini/GMOS 

1200 

20 

this work 

140423A 

3.258 

20.45 ±0.20 

> -1.44 

Fe 

Y 

Gemini/GMOS 

1200 

20 

this work 

140518A 

4.7055 

21.65 ±0.20 

> -1.06 

S 

Y 

Gemini/GMOS 

1200 

50 

this work 


Note. — List of the GRB-DLAs identified to date along GRB lines of sight. Missing metallicities are due to lack of metal line 
transitions or the low signal-to-noise of the spectra. When multiple measurements for the same line of sight are listed, the one 
derived from high resolution instrument is adopted in our analysis. Some values, either of A^hi oi" adopted from the 

literature reference listed in the last column. References: [IjGastro et al. (2003); [2]Vreeswijk et al. (2006); [3]Hjorth et al. (2003a); 
[4]Savaglio et al. (2002); [5]Fiore et al. (2005); [6]Shin et al. (2006); [7]Vreeswijk et al. (2004); [SjJakobsson et al. (2004); [9]Fynbo 
et al. (2009); [lOjLaskar et al. (2011); [11] Watson et al. (2006); [12]Berger et al. (2006); [13]D’Elia et al. (2007); [14]Prochaska 
et al. (2007); [15]Kawai et al. (2006); [16]Chen et al. (2007); [17]Fox et al. (2008); [ISjPenprase et al. (2006); [19]Chary et al. 
(2007); [20]Thone et al. (2010); [21]Ferrero et al. (2009); [22] Jakobsson et al. (2006); [23]Ruiz- Velasco et al. (2007); [24]Ledoux et al. 
(2009); [25]De Cia et al. (2011); [26]Prochaska et al. (2009); [27]Starling et al. (2009); [28]Page et al. (2009); [29] Patel et al. (2010); 
[30]D’Elia et al. (2011); [31]D’Avanzo et al. (2010); [32]Cenko et al. (2011); [33]Savaglio et al. (2012); [34]de Ugarte Postigo et al. 
(2012); [35]Skulad6ttir (2010); [36]Rau et al. (2010); [37]D’Elia et al. (2010); [38]Thone et al. (2013); [39]Cucchiara et al. (2011a); 
[40]Sparre et al. (2014); [41]D’Elia et al. (2014); [42]Kruhler et al. (2013); [43]Eriis et al. (2014); [44]Chornock et al. (2013). 

^Metallicities are relative to solar: [X/H\ = log[V/i7]^^^-log[V/i7]0. For the sub-DLA systems, with 10^® < Vhi < 

2x10^^ cm~^, we defer the reader to future paper. 

^On average over multiple spectra, see Appendix A 

^Value reported in the literature 
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GRB050820 - fflRES 
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Simulated Gemini 
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Relative Velocity (km s“^) 



-400 -200 0 200 400 


Fig. 1. — Left column shows weak metal lines identified in the Keck/HIRES {R = 30000) spectrum of GRB 050820A, which has 
true metallicity of Z/Zq = +0.17. Center and right columns show the same lines resampled at the typical X-Shooter {R ~ 6000 
in this case) and Gemini {R ~ 1200) resolutions respectively, assuming a S/N=10. Green line represent the continuum level. In 
case of metal rich systems, like this example, reliable abundances can still be derived from the X-Shooter data, though blending 
can be an issue, while only limits can be placed from the Gemini data. Lower resolution instruments are even more heavily 
affected, making it almost impossible to determine meaningful limits on the metal abundances. 
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GRB050730 - UVES 





-400 -200 0 


200 400 -400 -200 0 200 400 

Relative Velocity (km s“^) 


-400 -200 0 200 400 


Fig. 2 . — Same as Figure 1, this time for one of the lowest metallicity systems GRB 050730 (obtained with the VLT/UVES 
spectrograph, R ~ 40000) which has a true metallicity of ZJZq = +0.01. In this case it is obvious that the typical lines used 
for metal abundance measurements are also very narrow and almost disappear even at the X-Shooter resolution. This plot 
shows that measuring metallicity of the order of Z/Zq = +0.01 or less is very difficult and often only limits can be placed (e.g. 
S II1250 or Fe 11I6II). 
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Log A^hi 


Fig. 3 . — Column density distribution for our two DLA samples (QSO in red and GRB in blue). Vertical error bars are derived 
assuming poissonian distribution (95% confidence level) and we overplot in green (dashed line) the model by Noterdaeme et ah 
(2009). Despite the fact that the GRB-DLA sample is a factor of 4 smaller there is a clear overlap in the distribution around 
Vhi=21.5, while GRB-DLAs show a much larger number of dense systems compared with the QSO-DLA sample (see also 
Noterdaeme et ah 2014), indicating that GRBs are located, if not embedded, in very dense regions within their host galaxies, 
that can be the beacons of present star-formation. 


16 



Fig. 4. — Metallicity evolution with redshift for the GRB (red) and the QSO (grey) samples. Lower limits are indicated 
by upward triangles, while filled/open symbols indicate if these values come from high/low resolving power instruments. We 
perform a linear regression fit of the GRB-DLA data using the Schmitt survival analysis method, which keeps into account the 
censoring within the dataset (red dashed line). The shaded area represents the la error in the fitting parameters obtained using 
500 bootstrap iterations. A linear fit of the QSO-DLAs metallicity is marked by the dashed black line (see text for details). The 
GRB sample, despite the large scatter in [X/H], seems to probe an environment which slowly declines across the 2 : = 1.8 — 6 
redshift range. 
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Fig. 5 . — Adapted from Trenti et al. (2014). The colored lines represent two particular model for GRB hosts metallicity 
assuming a mild metallicity bias (p = 0.04, solid blue) and no bias at all (dotted green). Lines also indicate the trends for 
the top 95%, the median, and lower 20% of the simulated GRB host population. The mild bias model is broadly capable to 
reproduce the observed metallicity distribution of our biased sample of GRB-DLAs at a specific redshift, while the “no bias” is 
largely inconsistence with the data. 
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A. Description of relevant spectra 
A.l. GRB000926 

Our results confirm the finding published in the literature. The low-resolution data though allow to place 
only a lower limit based on the Zn ii lines (also consistent with Prochaska 2006). 

A.2. GRB011211 

Our examination of the FORS spectrum allows to place a lower limit based on the Si ii line of [X/H] > 
— 1.22, consistent with the result of (Prochaska 2006). 

A.3. GRB020124 

We derived similar Ahi value to the published work of Hjorth et al. (2003a), but the spectrum is too low 
S/N to provide an adequate metallicity measurement. 

A.4. GRB021004 

Extensive literature is present for this ORB sub-DLA (e.g. Fiore et al. 2005; Fynbo et al. 2005; Lazzati 
et al. 2006; Castro-Tirado et al. 2010). We report this GRB here for completeness. 

A.5. GRB030226 

We also notice the presence of an intervening system along this sightline at Zmt = 1.96236. Based on 
some isolated Iron lines we place a lower limit on this GRB-DLA consistent with the measurement of (Shin 
et al. 2006). Our iron-derived metallicity limit is consistent also with the measurement reported by (Schady 
et al. 2011), considering the large uncertainty in the HI estimate (0.3 dex). 

A.6. GRB030323 

Our analysis is consistent with the previous work by Vreeswijk et al. (2004) (including the different in 
the sulphur solar abundance between Asplund et al. (2009) and Grevesse & Sauval (1998)) and Prochaska 
(2006). 

A.7. GRB030429 

The spectra has been presented before Jakobsson et al. (2004). The lower limit in the metallicity is given 
by the few detected lines (e.g. Si ii). 

A.8. GRB050319 

The spectrum has been presented in Fynbo et al. (2009), we adopt the measurement from Laskar et al. 
(2011) as an upper limit based on Si ii 1526. 

A.9. GRB050401 

The spectrum was presented by Watson et al. (2006). We obtain similar results and adopt a value of 
[X/H] > —1.07 based on Zn ii line measurement. 

A.IO. GRB050505 

The metallicity we derived from sulfur measurements is similar to what has been presented by Berger 
et al. (2006), though we consider this to be a lower limit. 
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A.ll. GRB050730 


The high-resolution allows accurate determination of metallicity using different tracers. Our results pre- 
sented in Table 1 are consisting with the published values D’Elia et al. (2007) and Prochaska et al. (2008). 

A.12. GRB050820 

The measurement we obtained based on sulfur lines is consistent with the results from Jakobsson et al. 
(2006) and Prochaska et al. (2007). 

A.13. GRB050904 

As noted by Thone et al. (2013), the values reported in the literature (Kawai et al. 2006) is likely an 
overestimation due to blending of different lines, so we adopt their sulphur measurement, [S/H] > —1.0, as 
a lower limit for this GRB. 

A.14. GRB050908 

This is a sub-DLA system, so it is included for completeness. 

A.15. GRB050922C 

The values reported by Fox et al. (2008) and Prochaska et al. (2007) are consistent with our own mea- 
surement. 

A.16. GRB051111 

Only upper limit has been derived for this GRB. 

A.17. GRB060115 

These date are part of the compilation from Fynbo et al. (2009). Our estimated lower limit is mostly 
based on S ii lines. 

A.18. GRB060124 

This is a sub-DLA system, so it is included for completeness. 

A.19. GRB060206 

Our analysis is consistent with the results of Fynbo et al. (2006), but we consider our estimated value of 
[X/H] = —0.74 a lower limit, based on the strength of the same sulphur lines. 

A.20. GRB060210 

This GRB spectrum shows the presence of a nearby intervening systems at Zint = 3.817. Despite its 
presence a stringent lower limit of [X/H] > —0.83 can be placed based on iron and silicon lines. 

A.21. GRB060223 

The value from Ghary et al. (2007) has been adopted. 

A.22. GRB060510B 

We confirm the result from the literature based on sulfur lines Chary et al. (2007); Price et al. (2007). 


24 


A.23. GRB060522 

The value of -/Vhi from Chary et al. (2007) has been adopted, but there are no metal lines detected. 

A.24. GRB060526 

This is a sub-DLA system, so is included for completeness. 

A.25. GRB060605 

This is a sub-DLA system, so is included for completeness. 

A.26. GRB060607A 

This is a sub-DLA system, so is included for completeness. 

A.27. GRB060707 

We confirm the analysis of Jakobsson et al. (2006), but the low-resolution data allows only to place a 
conservative lower limit of [X/H] > —1.69 based on weak Fe ii transitions. 

A.28. GRB060714 

The spectrum has a higher S/N and similarly to GRB060707 we could derive a lower limit, but in this 
case we used Zn ii lines (Jakobsson et al. 2006). 

A.29. GRB060906 

This spectrum was presented by Jakobsson et al. (2006) and metallicity measurements were recently 
presented by Laskar et al. (2011). Using Sulfur lines we derive [X/H] > —1.72. 

A.30. GRB060926 

For this GRB we used likely saturated Zn ii lines and place a lower limit of [X/H] > —1.32. 

A.31. GRB060927 

This low S/N spectrum has been published before (Ruiz- Velasco et al. 2007). We derived more stringent 
limit then the ones derived by Laskar et al. (2011) using S ii saturated lines. 

A.32. GRB061110B 

This spectrum is part of the compilation of Fynbo et al. (2009). We obtained similar measurement than 
Laskar et al. (2011), but we consider this a lower limit. 

A.33. GRB070110 

This spectrum is part of the compilation of Fynbo et al. (2009), we performed our own metallicity 
measurement. 

A.34. GRB070411 

This is a sub-DLA system, so is included for completeness. 

A.35. GRB070506 

This spectrum is part of the compilation of Fynbo et al. (2009). We also identified an intervening system 
at z = 2.071. We derived a metallicity lower limit for the host from Zn ii lines. 
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A.36. GRB070721B 


This spectrum is part of the compilation of Fynbo et al. (2009). The low-resolution of this spectrum 
allows only a lower limit estimate from Si ii. 

A.37. GRB070802 

This spectrum is part of the compilation of Fynbo et al. (2009) and has been studied by Eliasdottir et al. 
(2009). The low-resolution of this spectrum allows only a lower limit estimate from blended Si ii: we are a 
bit more conservative than a previous work by Eliasdottir et al. (2009), but the two values are within 0.1 
dex, which is below the typical uncertainty for these low-resolution datasets. 

A.38. GRB071020 

This GRB is likely a sub-DLA system, so is included for completeness. 

A.39. GRB071031 

This spectrum was presented by Ledoux et al. (2009). Most of the useful metal lines are saturated (e.g. 
S II 1259). Using weak Fe ii features (e.g. Fe ii 1611) we obtained similar results [X/H] = —1.85 ±0.12 (see 
also Fox et al. 2008). 

AMO. GRB080210 

This GRB was observed with the FORS2 instrument at several different resolutions (see De Gia et al. 
2011). Our results, based on our analysis of the same FORS2 data shows similar results, though we consider 
these metallicity estimate a lower limit (see Section 3.1). 

A.41. GRB080310 

This GRB is likely a sub-DLA system, so is included for completeness. 

A.42. GRB080413A 

This spectrum was presented by Fynbo et al. (2009) and metallicity estimated by Ledoux et al. (2009). 
We derived only a lower limit based on zinc and nickel lines, which are also consistent with these previous 
works. 

A.43. GRB080607 

This GRB spectrum was presented by Prochaska et al. (2009) and it shows for the first time the clear 
presence in a GRB afterglow spectrum of H 2 molecular lines. Using different metallicity tracers, we retrieve 
similar values for the metallicity limits ([X/H] > —1.71). 

A.44. GRB080721 

This spectrum was presented by Starling et al. (2009). We derived a more conservative lower limit, using 
multiple a-elements lines. 

A.45. GRB080804 

The spectrum was presented by Fynbo et al. (2009). We determine the metallicity and relative depletion 
using Iron as well as non-refractory elements. 

A.46. GRB080810 

This GRB is likely a sub-DLA system, so is included for completeness. 
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A.47. GRB080913 


This GRB is likely a sub-DLA system, so is included for completeness. 

A.48. GRB081008 

This spectrum has been presented by D’Elia et al. (2011). We obtained similar metallicity estimate. We 
also analyzed a Gemini/GMOS spectrum from which we derived consistent lower limits. We used our UVES 
spectrum measurements. 

A.49. GRB090205 

This line of sight appears in D’Avanzo et al. (2010). In particular, we note that these author reported sulfur 
column abundance from blended profiles. We derived a similar metallicity lower limit of [X/H] > —0.57. 

A.50. GRB090323 

This GRB has a complex structure, with two absorbers within 600km s“^ Savaglio et al. (2012), one of 
which is a sub-DLA. Due to the complexity of this line of sight and the low-resolution of the VLT/FORS 
spectrum, we exclude GRB 090323 from our analysis. We report the metallicity in Table 1 for completeness. 

A.51. GRB090426 

This GRB is likely a sub-DLA system, so is included for completeness. 

A.52. GRB090516 

These data were presented by de Ugarte Postigo et al. (2012), but only equivalent widths were provided. 
Our analysis of the low-resolution VLT/FORS2 spectrum results in a lower limit of the metallicity of [X/H] > 

-1.36. 

A.53. GRB090809 

We measured a slightly higher metallicity than the one obtained by Skuladottir (2010): nevertheless our 
silicon-derived metallicity value (see Table 1) is well within 2cr from the one presented by these authors. 

A.54. GRB090812 

These data were presented by de Ugarte Postigo et al. (2012). We estimated a lower limit based on 
alpha-elements absorption features. 

A.55. GRB090926 

A VLT/FORS2 spectrum was presented by Rau et al. (2010), while an X-Shooter series of spectroscopic 
observations was presented by D’Elia et al. (2010). We obtained similar results, in particular the X-Shooter 
resolution allows stronger constrain on the metallicity. 

A.56. GRB100219 

This spectrum was presented by Thbne et al. (2013). We supplement this dataset with a Gemini/GMOS 
spectrum. We also reanalyzed the X-Shooter spectrum. Our measurement is slightly higher then the one 
presented by these authors, but within their error (0.2 dex). 

A.57. GRB100425 

We analyzed this dataset using the phase 3 products provided by the ESO database. Our metallicity 
lower limit is consistent with the one derived by Skuladottir (2010). 
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A.58. GRB110205 


This spectrum was already presented by Cucchiara et al. (2011a). 

A.59. GRB111008A 

This X-Shooter spectrum was presented by Sparre et al. (2014). Our Iron and sulfur metallicity estimate 
are consistent with the ones presented in this work. We report the former as metallicity tracer. 

A.60. GRB111107A 

This X-Shooter spectrum has not been presented in previous work. The signal-to-noise of the spectra is 
low, and therefore our metallicity measurement, based on saturated sulphur lines, has to be considered a 
lower limit. 

A.61. GRB120327A 

This spectrum has been published by D’Elia et al. (2014). Our sulphur measurement agrees with these 
authors. 

A.62. GRB120716A 

This spectrum has not been published in the literature before. The spectrum was obtained 2.6 days after 
the GRB has been discovered. Despite the low S/N we report a metallicity of [X/H] > —1.76 based on 
identified iron lines. 

A.63. GRB120815A 

This GRB spectrum has been published by Kriihler et al. (2013). Our analysis is consistent with these 
authors results. 

A.64. GRB120909A 

This X-Shooter spectrum has not been presented before. We derive metallicity measurements from weak 
Iron lines as well a-elements. 

A.65. GRB121024A 

This X-Shooter spectrum has been recently presente by Friis et al. (2014). We also identified multiple 
systems at Z\ = 2.3014 and an intervening one at Z 2 = 2.2977 (corresponding to 400 km s“^). Fine-structure 
transitions are identified in correspondence of Z 2 system, suggesting that the cloud at Zi is at large distance 
from the GRB radiation field. The broad, saturated Lya profiles makes hard to discern between the two 
components. Therefore we opted for considering them as one single absorber. The metallicity from zinc lines 
is reported in Table 1, while from iron lines we can infer a large depletion factor (see Friis et al. 2014, for a 
detailed analysis). 

A.66. GRB121201A 

This X-Shooter spectrum was preliminary presented in GCN only (Sanchez-Ramirez et al. 2012). This 
line of sight present a possible Lya emission line. Unfortunately the signal to noise is low (S/N< 3) and 
at the line redshift {z = 3.385) it is difficult to identify metal lines unambiguously. We therefore decide to 
report only the neutral hydrogen column density. 
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A.67. GRB130408A 


This GRB has been observed by our team Tanvir et al. (2013) as well as by VLT/X-Shooter Hjorth 
et al. (2013). We analyzed both spectra and report the metallicity measurements from both analysis. The 
X-shooter derived value is plotted in Figured. 

A.68. GRB130505A 

We present our Gemini spectrum for the first time. We identified several absorption lines, and we were 
able to determine a lower limit on the metallicity based on iron weak transition (Fe 11I6O8). 

A.69. GRB130606A 

This spectrum has been published by Ghornock et al. (2013). The line of sight present signature of a 
likely sub-DLA system, therefore we present this value here for completeness. 

A.70. GRB140226A 

This spectrum was obtained by the Keck/LRIS instrument. Several metal lines have been identified, but 
the low-resolution of this instrument allows to place only a metallicity lower limit from Iron and Sulphur 
lines. 

A.71. GRB140311A 

This spectrum was obtained by our collaboration and present few metal absorption features. We were only 
able to place a lower limit on the metallicity based on nickel lines, since sulphur lines seems contaminated 
by other lines. 

A.72. GRB140419A 

This is a sub-DLA system, so is included for completeness. 

A.73. GRB140423A 

We present our analysis on our Gemini/GMOS spectrum (Tanvir et al. 2014). Based on a-elements lines 
and weak iron lines (e.g. Fe ii 1608) we were able to place a lower limit on the metallicity. 

A.74. GRB140518A 

This data have not been published before, but preliminary analysis appears in Ghornock et al. (2014). 
Our metallicity limit comes from saturated sulphur transitions. 
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